Spin Filtered Edge States and Quantum Hall Effect in Graphene 
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Electron edge states in graphene in the Quantum Hall effect regime can carry both charge and 
spin. We show that spin splitting of the zeroth Landau level gives rise to counterpropagating modes 
with opposite spin polarization. These chiral spin modes lead to a rich variety of spin current states, 
depending on the spin flip rate. A method to control the latter locally is proposed. We estimate 
Zeeman spin splitting enhanced by exchange, and obtain a spin gap of a few hundred Kelvin. 



A new electron system with low carrier density and 
high mobility was recently realized in two-dimensional 
graphene 1] . By varying the carrier density with a gate 
one can explore a range of interesting state^ in particu- 
lar the anomalous quantum Hall effect '3,'^ (QHE). In 
contrast to the well-known integer QHE in silicon MOS- 
FETs 3 the QHE in graphene occurs at half-integer mul- 
tiples of 4, the degeneracy due to spin and orbit. This 
has been called the half-integer QHE. The unusually large 
Landau level spacing makes QHE in graphene observable 
at temperatures of 100 K and higher. 

Here we explore the spin effects in graphene QHE. In 
the presence of Zeeman splitting transport in graphene is 
described by an unusual set of edge states which we shall 
call chiral spin edge states. These states are reminiscent 
of the ordinary QHE edge states @ , but can propagate in 
opposite directions for opposite spin polarizations. (As 
shown in 0, similar states can arise due to spin-orbit 
coupling in the absence of magnetic field. However, the 
weakness of spin-orbital effects makes the corresponding 
spin gap quite small.) The chiral spin edge modes can 
be used to realize an interesting spin transport regime, 
in which spin and charge currents can be controlled inde- 
pendently. Observation of these phenomena is facilitated 
by fairly large magnitude of the spin gap. The gap is en- 
hanced due to electron correlation and exchange, and can 
reach a few hundred Kelvin for realistic magnetic field. 

The half-integer QHE in graphene was interpreted 
in terms of a quantum anomaly of the zeroth Landau 
level Alternatively, these properties are easily un- 
derstood from the edge states viewpoint, similar to the 
usual QHE. This was done is Ref. (s| using numerical 
treatment of the zigzag edge. Here we present a contin- 
uum description of the edge states, using the massless 
Dirac model which provides a good approximation for 
Carbon 7r-electron band near its center. We reduce the 
problem to a one-dimensional Schrodinger equation with 
a potential which depends on the boundary type. By 
comparing the behavior for armchair and zigzag bound- 
ary, we show that the energy spectrum properties near 
the edge are universal and imply the half-integer QHE. 

To interpret the half-integer QHE, let us inspect the 
energies of the first few Landau levels (LL) obtained for 
an armchair boundary (Fig. ^a)). First we ignore elec- 
tron spin. In the bulk the LL's are doubly degenerate. 
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FIG. 1: (a) Graphene energy spectrum near the armchair 
boundary obtained from Dirac model, Eq.Q. The boundary 
condition, Eq.l|^, lifts the K, K' degeneracy. The odd integer 
numbers of edge modes lead to the half-integer QHE. (b) Spin- 
split graphene edge states: the blue (red) curves represent 
the spin up (spin down) states. These states propagate in 
opposite directions at zero energy. 

due to two species of Dirac particles located near K and 
K' , the inequivalent corners of the first Brillouin zone. 
We note that the zeroth LL splits into two levels with 
positive and negative energies. In contrast, the behavior 
of the edge states associated with higher LL's is more con- 
ventional 01 : the energies of positive (negative) LL's in- 
crease (decrease) as one approaches the boundary. Hence 
in the spinless case the number of edge states can take 
only odd integer values and the Hall conductivity is odd 
integer in units of e^/ft.. For example, when the chemi- 
cal potential is between the n = —1 and n — —2 LL's, 
there are three branches of active edge states: two of 
them derived from the LL with n = — 1 and one from 
the LL with n = 0. As a result, although each Landau 
level filling factor is an integer, the conductance at QHE 
plateaus is half-integer in units of Ae^/h which accounts 
for both the K, K' and spin degeneracy. 

This behavior is modified in an interesting way by the 
spin splitting of LL's (Fig. ^b)). When the chemical 



potential /i lies in the interval 



< /i < iAs, the 



zeroth LL is spin polarized, with only spin down states 
being filled. However, there exists a branch of up-spin 
edge states going to negative energies. The states of this 
branch with e < n will be filled and will contribute to 
transport on equal footing with the down-spin states. 
Notably, the up-spin and down-spin states have oppo- 
site chiralities, i.e. they propagate in opposite directions. 
These states carry opposite charge currents but equal 
spin currents. As a result, the edge transport in the spin 
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gap, — < fi < ^As, can be spin filtered. 

The chiral spin edge states found here are similar to 
those predicted by Kane and Mele 0, who considered 
spin-orbit coupling in graphene in the absence of Zeeman 
field. They obtain a spin gap which is solely due to the 
spin-orbit coupling and is estimated as Aso ~ IK 
The small gap poses strong constraints on the amount 
of disorder and on temperature, hindering experimental 
tests of the interesting predictions of Ref. p. In con- 
trast, the magnetic field-induced gap is quite large: the 
Zeeman energy of Aj = gusB ~ 15 K at B ~ 10 T is fur- 
ther enhanced by exchange interaction, putting it in an 
experimentally convenient range of hundreds of Kelvin. 
Furthermore, one can alter the character of spin flip scat- 
tering at the edge to induce backscattering among the 
spin polarized edge modes, enabling experimental obser- 
vation of a number of novel transport effects (see below). 

To analyze the edge states, we employ the Dirac 
model 01 which describes the low-lying states as linear 
combinations of four zero-energy Bloch functions with 
slow varying envelope functions uk, vk, uk', vk'- Here 
u and V correspond to the wave function components on 
two inequivalent atomic sites A, B, and the subscripts 
K , K' denote the inequivalent Dirac points. The effec- 
tive low-energy Hamiltonian, written near each of the 
points if, K' in terms of uk,k' , vk,k' , is of the form 
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where p± = px±ipy, Pf^ =Pfj_-^Ai^, andvo « 8xl0''cm/s 
is Fermi velocity. The energy spectrum in the bulk is 
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sgn(n)|n|i/2g^^ = hvo {2eB /hc)^^^ (2) 



with integer n. For typical magnetic field of 10 T, the 
lowest LL separation is estimated to be quite large. 
So ^ El — Eq ~ 1000 K. We now consider a graphene 
sample with an armchair edge parallel to the x axis, us- 
ing Landau gauge = —By, Ay = 0. We eliminate 
the V components of the wave function and consider the 
eigenvalue equations for the u components: 

I {pI + {v- Vo? + 1) - {E/e^fuK, (3) 
\ {pI + {y~ yo? ~ 1) UK' = {E/eo)^UK', (4) 
where yo = —Px- Here x, y and px, Py are measured in 

1 /2 

the units of is — {Tic/eB) and U/iB- 

To obtain the energy spectrum near the edge, one 
needs to supplement Eqs.lO,l@J with suitable boundary 
conditions. We obtain the latter from the tight-binding 
model assuming that it is valid up to the very last row 
near the boundary. In the armchair case the boundary 
condition requires the wave function to vanish at both A 
and B lattice sites along the line y = 0. For the envelope 



functions u, v, taken at j/ = 0, this means 

uk^uk', vk^vk' (5) 

(see Ref. 01 for a more general analysis). 

Instead of the problem in a halfplane with boundary 
conditions |(SJ), it is more convenient to consider Eq.© 
on the negative half-axis y < 0, and Eq.Q on the pos- 
itive half-axis y > Q. Then the first boundary condi- 
tion in Eq.© implies continuity of the wavefunction at 
?/ = 0, while the second condition implies continuity of 
the derivative dujdy. Therefore the problem reduces to 
calculating energy levels in the potential 



^(y) = ^(l2/|-yo)'-^sgn(y). 



(6) 



defined on the entire line. Numerical solution of this 
problem gives the energy dispersion shown in Fig. ^ 

We now briefly discuss the zigzag edge, another very 
common graphene edge type. In the absence of magnetic 
field this boundary supports a band of low-energy sur- 
face states [13, El • Wavefunctions of these states decay 
away from the boundary on a length scale of the order of 
lattice spacing, which is typically much shorter than the 
magnetic length Ib- Therefore, the surface states should 
not be sensitive to the magnetic field. This conclusion is 
supported by the numerical study f^] . 

In contrast to the surface states, special for zigzag edge, 
the LL behavior near the edge is universal and depends 
on the boundary type only weakly. In the Dirac equation 
framework, the boundary condition for a zigzag edge is 
that either A ox B components of the envelope function 
vanish at the boundary, depending on the orientation of 
the edge. This yields two possible boundary conditions: 
either uk = uk> = or vk = vk' = at the edge. Find- 
ing energy spectrum in this case requires solving Landau 
level problem with a hard-wall boundary condition, fa- 
miliar from the usual QHE j3] . The energy spectrum for 
the zigzag edge, given by the square root of the disper- 
sion curves found in Ref. is qualitatively similar to 
the one in Fig. These results are consistent with the 
numerical analysis of Ref. 0^ . 

Turning to the analysis of the spin gap, let us consider 
the exchange energy for half-filled spin-degenerate Lan- 
dau level. In the Hartree-Fock approximation, it is given 
by spin-dependent correlation energy: 
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V{r -r')gij{r -r')<fr<rr' (7) 



with gij{r — r') = {ni{r)nj{r')) the pair correlation fimc- 
tions. We consider Coulomb interaction of the form 
Vir) — / nr with the dielectric constant [3 describ- 
ing screening due to the filled states, e < 0. 

To describe correlations, we use the result for f ully 
filled Landau level [13, 9o{r) = 1 - e-'^'/^^^, setting [ig 

eB 

g++{r) = nlrilgoir), g^^{r) = n?_n^go(0> = 
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where n± are the occupation fractions of the up and 
down spin states, n+ + n_ = 1. We model the Coulomb 
interaction-induced correlations of particles with oppo- 
site spin, which are absent for noninteracting electrons, 

by (7-1 = n-|-n_ng(l — ae"'' /^^^), where a < 1 describes 

relative strength of Coulomb and exchange correlations. 
This approach yields 



^-An,n., A={^y^'^il-a). (8) 



exchange 



To estimate spin polarization, we consider Gibbs free en- 
ergy per particle for the half-filled spin degenerate Lan- 
dau level, G/N = E-TS, where 

S = — T rii Inn-i, E ~ —Ez{n-\- — n_) — ^n-|_n_, 

with Ez the Zeeman energy. The ground state, deter- 
mined by 5G = with + = 1, satisfies 



rin(n+/n_) - Ez + A{n^ - n+) = 0. 



(9) 



In the absence of exchange interaction, A — 0, this would 
give n^/n+ = e^^^l^ . The more realistic case of ex- 
change A large compared to Ez can be analyzed by set- 
ting Ez = 0. In this case we have a phase transition at 
Tc = \A. At low temperature, T ^ Tc, the concentra- 
tions n± satisfy n-/nj^ = e.~^l^ ^ i.e. the spin gap at 
T can be estimated as A = A. 

To obtain numerical value of A, we use the RPA esti- 
mate of the screening function 14], with /Tiv « 2.7: 



K = 1 + 27re^ 



1 

Afiv 



5.24. 



Comparing to the LL separation £q, Eq.(|2Jl, 
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■(1 -a)£o« 0.456 •(l-a)eo- 



(10) 



(11) 



Note that A is proportional to the Dirac LL energy Eq 
and scales as B^/^, in contrast to the Zeeman energy. If 
we use a = 0, i.e. ignore correlations of electrons with 
opposite spins, we obtain A ~ 450 K for B ~ 10 T. (This 
approximation, while giving correct order of magnitude, 
may somewhat overestimate the exchange contribution.) 

We now discuss possible experimental tests of the chi- 
ral spin edge states, using the four-terminal device shown 
in Figs. I2I3I We assume that each contact injects both 
spin polarizations with the same voltage 14, i.e. full spin 
mixing takes place in the leads. The charge current flow- 
ing out of the fc-th contact is given by 



4' = E5fcfc' iVk-Vk-), 



(12) 



where g^k' is the Landauer conductance of the edge chan- 
nel connecting contacts k and fc', and Vk are voltages on 
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FIG. 2: Four-terminal device with chiral spin edge states (no 
backscattering) . In response to charge current of I'^ = 2iv = 
4:e'^V/h between reservoirs 1 and 2, pure spin current of P = 
2iv flows between 4 and 3, while Hall voltage is zero. 
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FIG. 3: Same as in Fig. |5| with backscattering between 
edge states induced by altering local spin-flip rate. Strong 
backscattering gives rise to asymmetric current flow with the 
current between 1 and 2 fully spin polarized. Inset shows how 
Hall voltage probe can be used to detect spin current. 

the contacts. In the presence of backscattering among 
the edge states, the conductance gkk' can be expressed 
in terms of transmission matrix Tkk' of the edge chan- 
nels 01, with gkk' = Tkk'e^/h. The voltage probes are 
defined by the condition = and Eq. p2ll is solved for 
the current and voltage at each contact 13]. 

The spin current is then determined as follows. The 
spin current flowing from the reservoir k to reser- 
voir k' is given by the sum of outgoing current of 
e'^Vk/h, reflected current (1 — Tkk')e^Vk/h and the cur- 
rent from fc'-th reservoir, Tkk'c'^Vk' /h. We obtain I^^,, = 
±[{2 — Tkk')Vk +Tkk'Vk']e'^/h (the three contributions 
to the spin current are of the same sign, plus if the chan- 
nel fc — > A:' is spin up and minus if it is spin down) . Thus 
the total spin current flowing out of the A:-th contact is 



kk' 



k' 



£kk'gkk'{Vk - Vk'), 



(13) 



where ekk' = —£k'k equals -1-1 (—1) when the current 
from k to k' is carried by spin up (spin down) electrons. 

The general relations (|12l) , H13|l become more transpar- 
ent in two simple limiting cases, illustrated in Figs.|21ISl 
First we inspect the clean limit with no backscattering 
at the edge (Fig. Applying voltage of 21^ between 
contacts 1 and 2 we obtain pure spin current across the 
sample equal to 2iy, where iy = e^Vjh. There is also 
charge current of 2iy between 1 and 2 and no Hall volt- 
age, yielding p^^ = hj^e^. 

Alternatively, let us introduce strong backscattering in 
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the shoulder connecting 1 and 4, such that the resistance 
between 1 and 4 is infinite, while all the other shoulders 
are clean (Fig. 13). Such a system acts as a "spin filter". 
Indeed, the current flowing from 1 to 2 is spin polarized. 
The Hall voltage in this case is nonzero, equal to the half 
of the voltage between 1 and 2. 

Interestingly, the Hall voltage measures spin current 
rather than charge current. This is illustrated by Fig.|3| 
inset which depicts currents zi, Z2 that flow into the volt- 
age probes, equilibrate and flow out {i'l 2 — 5*1,2)- The 
Hall voltage is proportional to the incoming spin current: 

Vi-V2 = ^{ii - i2)h/e^ = -Is ■ h/2e\ (14) 

Such voltage probes can be attached to the top or bottom 
leads in Fig.|21l2|to verify the presence of spin currents. 
(We caution that currents ii, 12 are affected by the probes 
and should be calculated self-consistently using ^^.) 

We now discuss spin-flip scattering mechanisms. One 
can make qualitative observations which (a) show that 
spin-flip scattering is strongly suppressed and (b) sug- 
gest a way to manipulate it. We assume that the main 
source of the spin-flip scattering is spin-orbit (SO) inter- 
action, which includes the intrinsic term proportional to 
Gz and the Rashba term 0. The former interaction is 
ineffective when electron spins are polarized along the z 
axis, amounting merely to an energy shift of the order of 
^SO ~ 1 K, without the up-down spin mixing. At the 
same time, the Rashba term turns out to be extremely 
small, \ji ^ 0.5 mK making spin-flips negligible. 

The situation changes, however, when an in-plane 
magnetic field is applied. Now the intrinsic SO inter- 
action admixes the spin projections along the tilted field, 
leading to a small avoided crossing of the chiral spin edge 
states of magnitude A50. Using a local gate one can de- 
tune from this crossing by Ai? ^ A50. In this case, the 
spin-flip scattering probability is estimated as Al^/Ai?^, 
which gives a factor of 10~^ for IS.E ~ 10 K. The spin-flip 
scattering rate will be further reduced if the disorder at 
the edge has a short length scale compared to ,^_b (recent 
STM experiments 0,^3 report a ~ Inm). For typical 
magnetic length, Ib ~ 10 nm at B ~ 10 T, this gives spin- 
flip scattering suppressed by a factor of a'^/lg ~ 10~^. 
These observations show that backscattering at the edge 
should be strongly suppressed and suggest that by gat- 
ing the QHE edge in the presence of in-plane magnetic 
field one can control local spin-fiip rates independently. 
In particular, by tuning AE to zero in one shoulder of 
the four terminal device, one could realize the situation 
depicted in Fig. 13 

Up to now we have focused on the LL spin splitting 
and ignored the orbital K, K' degeneracy. However, the 
exchange effects must lift the K, K' degeneracy as well. 
Moreover, the K, K' exchange gap should be of the same 
size as the spin gap A since electron interaction has the 
same strength for different spin and valley species {SU (4) 



symmetry). The spontaneously appearing order param- 
eter will have orientation fixed by weak nonsymmetric 
interactions, the difference between the intra and inter- 
sublattice coupling for K, K', and the Zeeman energy for 
spin splitting. The spin filter effects discussed above as- 
sume that spin splitting occurs before orbital splitting. 

We conclude by addressing the experimental situa- 
tion. Our theory predicts a plateau of Rxy = be- 
tween eV = ±^A which was observed in a recent ex- 
periment 0. This paper also reports full lifting of the 
four-fould degeneracy of the zeroth LL at high magnetic 
fields B > lOT which is consistent with the predicted 
strength of K, K' exchange. We attribute the suppres- 
sion of the zeroth LL splitting at lower magnetic fields 
to the disorder present in the currently available sam- 
ples |23|. The observed value of Rxx &t B = 25T is 
about 40kri. According to our picture, this means that 
the chiral edge states are not localized and transmission 
coefficients of the edge channels are about 1/3. Thus a 
spin Hall current may already be present in this case. 

After the completion of this work, we became aware of 
Ref. [U where a continuum description of the edge states 
given in the first part of this paper is also reported. 
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